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Most hadronic particles are resonances: for example, the p meson appears as a resonance in the elastic scattering 
of two pions. A method by Liischer enables one to measure the properties of the resonance particles from finite 
lattices. We present a more general method which includes scattering processes where the total momentum of 
the particles is non-zero. The main advantage is that the resonance scattering can be observed in a considerably 
smaller spatial volume. We test the method with a simple 3+1 dimensional spin model, and find excellent 
agreement between the zero momentum and the non-zero momentum scattering sectors. 



1. INTRODUCTION 

Lattice simulations have been very succesful in 
resolving many low energy properties of hadronic 
particles. However, many low energy hadrons are 
resonances, which cannot be fully described with- 
out scattering into asymptotic states. The best 
known example is the p meson which appears in 
the elastic 7T7t L = 1, I = 1 scattering channel. 

Elastic scattering involves at least four asymp- 
totic states, which, in the case of QCD, are them- 
selves bound states of fundamental quark and 
gluon fields. Because the lattice simulations are 
restricted to small (periodic) physical volumes, 
the direct observation of free asymptotic states 
is not possible. At best, the lattice size is a few 
times the interaction length. 

This problem has been addressed in a series of 
papers by Liischer [1], giving a relation between 
the infinite volume elastic scattering phase shift 
and the 2-particle energy spectrum measured on 
finite periodic lattices. The finite volume is used 
to probe the system, and simulations with several 
volumes are needed to obtain full knowledge of 
the phase shift function. This method has been 
applied to various test models in 2-4 dimensions 
[2-5]. 

Liischer's method requires that the scattering 
particles have zero total momentum P = 0. We 
have generalized the method to arbitary total mo- 
mentum P [6]. This method has several desirable 
features: 

• When P ^ the resonance can be observed 
in smaller volumes than with P = 0. 



• P = and ^ sectors can be measured si- 
multaneously; depending on the model, of- 
ten with negligible cost in computer time. 

• For a given lattice size, the P = and ^ 
methods yield the phase shift 6i(p) with a 
different argument p - two data points from 
a single run. 

To illustrate the difference in the required lat- 
tice size, consider a 1-dimensional system contain- 
ing a light particle <f>, and a heavy resonance p. 
Let us first consider the energy levels in the non- 
interactmg case. With lattice size L, the lattice 
momenta are quantized to p = n2Tr/L, n £ Z. 
When we later turn on a small p<f><f> interaction, 
the resonance appears when the energies of the p 
and 2<f> -states are equal: W p = For con- 

creteness, let us also choose m p = |m^. This 
choice is by no means unique, any value satisfy- 
ing the limits 2m^ < m p < 4ra^ would do as 
well. 

(a) P = sector: the ^-particles have opposite 
momenta, the smallest value is p^ = (2tt/ L), and 
P P = 0- 




Using the "would-be resonance" condition W p = 
W2(j>, we see that the required lattice size is L ps 
10.5/m^. 

(b) P zfz 0: now we can choose p^ i = (2tt / L) and 
p^ 2 = 0, so that p p = (2-7T / L) and 
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This gives L 



4.5/m^ at the resonance. The 



ratio of volumes in cases (a) and (b) is y — 

(ij^) 3 Rj 12.7! The non-zero momentum sector 
gains more than an order of magnitude over the 
zero momentum case. 
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Even in _P ^ case it is convenient to work in 
the the center of mass frame of the scattering 
particles. In this frame the 2-particle energy is, 
in terms of the lattice (laboratory) frame energy 



Won, 2^ - 
and the momenta of the particles 
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In fig. 1, we show the 2-particle center of mass 
frame energy spectrum W'i$ both for P = and 
P = 2-7T / L sectors. The lattice sizes calculated 
in the previous example are indicated by arrows. 
When we turn the (small) p<f)(f> interaction on, the 
points where the 2<f> levels cross the p levels turn 
into avoided level crossings , and p becomes un- 
stable. 

2. PHASE SHIFT 

The phase shift 6j parametrizes the partial 
wave decomposition of the scattering amplitude: 
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In 3 + 1 dimensions the connection between the 2- 
particle energy levels and S is quite involved, and 
we present only the final formula here [6]. For 
concreteness, let P = 2ttx/L. The energy W C m,2^ 
is measured, and the CM-frame momentum p CTn 
is calculated from 
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Figure 1. The free particle center of mass energy 
levels for the P = (dashed lines) and P = 2tt / L 
(solid lines) sectors. 
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In the s-wave (I = 0) sector, the relation between 
p CTn and the phase shift 8$ is given by 
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where $ is a continuous function defined by 

tM-m) = J q f\, $(o) = o, (6) 



z 00 (i;q 2 ) 

with 7 = (1 - v' 2 )- 1 !' 2 = (1 + P' 2 /W' 2 m ) 1 l' 2 . Func- 
tion Zoo is a generalized zeta function, given by 
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where the set A = {r G R 3 |ri=7 _1 (rii + 
1/2), r2=«2, f*3=«3; n G Z 3 }. The expansion of 
Zoo is convergent when Res > 3/2, but it can be 
analytically continued to s = 1 [6]. When P — ► 0, 
7 — v 1 and A — ► Z 3 , and Liischer's original for- 
mula [1] is recovered. 

3. SIMULATIONS AND RESULTS 

We test the method with a simple 3+1 dimen- 
sional spin model containing two coupled Ising 
spins 
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Figure 2. The center of mass energy levels for sectors P = (top row) and P = 2tt / L (bottom) for cases 
A, B and C (see table 1). 



S = ^2 [-K,p<f> x <f>y - K p p x p y + gp x <f> x <f>y] , (8) 

where <f>,p = ±1. Choosing the coupling con- 
stants suitably p appears as a resonance in s-wave 
<f><f> scattering. Let us now define operators 

0(P,p;t) = 4>(P + p;t)4>(-p;t), (9) 

where <f> is the spatial Fourier transform of <f>. We 
calculate the phase shift as follows: (i) we mea- 
sure the cross-correlation matrices of operators O 
for both P = and P = 2tt / L separately; (ii) the 
exponential falloff of the eigenvalues of the matri- 
ces yields the energy levels; (iii) using eqs. (4-7) 
we obtain 6o(p C m) [6]. 

We use 3 sets of couplings (A, B, C), given 
in table 1. In fig. 2 we show the center of mass 
energy levels for both of the P sectors. In case A 
the coupling g = and p is stable; in cases B and 



C g ^ and the crossings of the p and 2<f> levels 
turn into avoided level crossings. In the P = 
sector even the lowest 2<f> energy level increases 
when L decreases. This is due to the repulsive 
\(j> 4 interaction in the scaling limit of the Ising 
model, and in cases A and B it causes the level 

Table 1 



The couplings of runs A, B and C, and the mea- 
sured scattering parameters. 
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0.0742 
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0.07075 
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0.0718 


0.0665 
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Figure 3. The phase shifts for cases A (g = 0), B 
(0.008) and C (0.021). Filled and open symbols 
correspond to P = and 2tt/L, respectively. 



crossing at L ps 8. However, these points are 
plaqued by severe finite volume effects, and the 
first useful crossing occurs only around L = 37. 
This shifts down to L ps 16 when P ^ 0. 

In fig. 3 we show the phase shifts for P = 
and 2tt/L for A, B and C. The repulsive A^ 4 in- 
teraction is evidenced by the negative slope (oc 
scattering length) of 5 near p = 0. In B and C 
the resonance appears as a rapid increase of 6 



near W c] 



The width of the resonance is 



excellent consistency of P = and 2tt / L sectors 
(filled and open symbols). 

The continuous lines in fig. 3 are perturbative 
fits to the So data: 
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The fitted parameters Xr, gR and m p are listed in 
table 1. Inverting eqs. (4-7) we obtain the energy 
levels corresponding to (5 pert , as shown in fig. 2 
with continuous lines. Note that all the levels 
are fitted simultaneously. Finally, the resonance 
width is calculated from 
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Conclusions 

We have measured resonance scattering param- 
eters to a high accuracy using Monte Carlo meth- 
ods in 3+1 dimensions. The sectors P = and 
2-7T / L are completely compatible within statisti- 
cal errors; however, level crossings occur in the 
P = 2-7T / L sector in much smaller volumes. The 
sectors complement each other: for a fixed vol- 
ume, they give S(p) with different argument p; 
and both sectors can be calculated simultane- 
ously, often at little extra cost. 
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proportional to the slope of 6 at 7r/2. Note the 



